Photon-assisted Landau-Zener transition: Role of coherent superposition states 
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We investigate a Landau-Zener (LZ) transition process modeled by a quantum two-level system 
(TLS) coupled to a photon mode when the bias energy is varied linearly in time. The initial 
state of the photon field is assumed to be a superposition of coherent states, leading to a more 
intricate LZ transition. Applying the rotating- wave approximation (RWA), analytical results are 
obtained revealing the enhancement of the LZ probability by increasing the average photon number. 
We also consider the creation of entanglement and the change of photon statistics during the LZ 
process. Without the RWA, we find some qualitative differences of the LZ dynamics from the RWA 
results, e.g., the average photon number no longer monotonically enhances the LZ probability. The 
ramifications and implications of these results are explored. 
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I. INTRODUCTION 

Landau-Zener (LZ) transitions involve a quantum two- 
level system (TLS) with a constant coupling strength A 
between two adiabatic energy levels. A control parame- 
ter is swept at a constant velocity v, so that an avoided 
crossing of energy levels occurs, and provides the prob- 
ability that the system will stay in an adiabatic state. 
LZ transitions have attracted considerable attention the- 
oretically (see, e.g., Refs. [1-6]) and experimentally (see, 
e.g., Refs. [7-9]). 

In a variety of physical areas, LZ processes play an 
important role, e.g., in artificial atoms [10] and Bose- 
Einstein condensates in optical lattices [11]. Especially 
in superconducting circuits which can behave like con- 
trollable quantum TLSs [12-15], LZ and Landau-Zener- 
Stiickelberg (LZS) problems have been studied by several 
groups [1, 16-21]. The standard LZ problem for an iso- 
lated TLS can be solved exactly. For some many-level 
systems, the LZ transition probability can also be cal- 
culated exactly for some initial states [22-24] . LZ tran- 
sitions controlled by classical fields are considered in 
Ref. [25], and in a quantum photon field the authors of 
Ref. [26] found that varying the LZ sweep rate produces 
collapses and revivals of the coherent field amplitude. 

Naturally, a quantum TLS is influenced by its environ- 
ment, and therefore there have been many studies about 
the dissipative LZ problem. Exact results are available at 
zero temperature [27, 28] , and various numerical methods 
have been employed to study the cases at finite temper- 
atures [29-31]. The nonmonotonic dependence of the LZ 
probability on the sweep velocity was studied in Ref. [29] 
using numerical methods. Environment parameters, such 
as temperature, can exponentially enhance the coherent 
oscillations generated at a LZ transition [31]. 
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It is interesting to replace the classical coupling A by a 
fully quantum-field coupling; then the TLS and the field 
form a whole composite quantum system. In this paper 
we consider a quantum TLS coupled to a photon mode 
when the bias energy is varied linearly in time. 

Coherent superpositions of coherent states, like 
\ip (0)) h = (| a) + e w |— a)) /N$, with the normalization 

constant Ng = 2(1 + cos#e~ 2 l Q ), have attracted ex- 
tensive interest as a distinct class of nonclassical states 
with interesting properties. For a large amplitude a, 
these can be interpreted as quantum superpositions of 
two macroscopically-distinguishabie states, the so-called 
Schrodinger cat states. Such states can be prepared in 
various systems and play an important role in funda- 
mental tests of quantum theory and in many quantum- 
information-processing tasks [32-34] , including quantum 
computation [35] , quantum teleportation [36] , and preci- 
sion measurements [37, 38] . 

We aim to discover the effect of the initial superposi- 
tion of coherent states on the LZ transition. The increas- 
ing average photon number may enhance the LZ proba- 
bility. We also focus on the effect of the LZ process on 
the quantum properties of the whole system, including 
entanglement creation and changing the photon distri- 
bution. 

By applying a rotating- wave approximation (RWA) , we 
obtain analytical results which reveal the enhancement 
of the LZ probability when increasing the average pho- 
ton number. Whereas, without the RWA we find some 
qualitative differences of the LZ dynamics from the RWA 
results; e.g., there are two stages of the LZ transition and 
the final LZ probability no longer monotonically depends 
on the average photon number. 

This paper is organized as follows. In Sec. II, we in- 
troduce the standard LZ model and the quantized LZ 
model considered in this paper. In Sec. Ill, by employing 
the RWA, we analytically calculate the LZ probability, 
the entanglement between the TLS and the field, and the 
photon statistics characterized by the Mandel parameter 
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Q. Numerical results are shown in this section to confirm 
the analytical results. In Sec. IV, without the RWA, nu- 
merical and analytical results are given to compare with 
the RWA results. The thermal state of the photon field 
is also considered, in order to compare with the case of 
a superposition of coherent states in which the photon 
distribution is Poissonian. Finally, we present the con- 
clusions. 



II. HAMILTONIAN OF THE LANDAU-ZENER 
TRANSITION 

Let us first briefly gather together several results that 
will be used in this work. The standard LZ problem, for 
an isolated quantum TLS driven by a classical field, is 
described by the Hamiltonian 



Hamiltonian becomes 



H 



vt 



(1) 



in terms of the Pauli matrices a XyZ , and a x = a + + cr_ 
(h = 1 is assumed throughout). Let the states |f) and |t) 
denote the eigenstates of a z , i.e., the so-called diabatic 
states with energies ±vt/2 which cross at t = 0. The 
constant v is the sweep velocity, by which the energies 
of the diabatic states cross. The coupling A denotes the 
interaction between the two diabatic states, which are 
chosen to be positive and time independent. For A^O, 
the diabatic states are not eigenstates of the Hamilto- 
nian in Eq. (1), and the avoided- level crossing appears be- 
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tween the adiabatic energies E± 

at t = 0. Thus, generally, a population transfer is in- 
duced. Asymptotically, for times \t\ S> A/v, the di- 
abatic states coincide with the adiabatic states. The 
LZ problem asks for the probability of the TLS end- 
ing up in the initially unoccupied level, and is given by 
Po.lz = 1 — cxp(— 7rA 2 /2w), which is an exact result for 
all A and v. In the adiabatic limit A 2 /v ^> 1, i.e., when 
the sweep occurs slowly enough, Po.lz will saturate at 
1, which implies that the transfer of population between 
the adiabatic eigenstates is prevented by the splitting A. 

In this work we consider a quantized LZ Hamiltonian 
describing the coupling of a quantum TLS to a single- 
photon-field mode. The Hamiltonian reads 
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where the operator a (cv) annihilates (creates) a photon 
in the field mode with frequency w, and the energy bias of 
the TLS is denoted by loq. In this paper, the resonance 
case wo = u> is considered. Hence we rewrite the total 
Hamiltonian in a rotating frame defined by the operator 
N = a^a + o~ z /2, at the frequency uj. 

In the weak-coupling regime, where the coupling is at 
least an order of magnitude less than the energy fre- 
quency, i.e., A < O.lw, one can employ a RWA, and the 
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Now the system is modeled in terms of a time-dependent 
Jaynes-Cummings Hamiltonian [26, 39, 40], which can 
also be called the Landau-Zener- Jaynes-Cummings (LZ- 
JC) model. Note that the operator N = a^a + a z /2 is 
conserved by this Hamiltonian. 

When considering a finite detuning 8ui = ujq — uj, the 
RWA is usually justified in the condition that \6lj\ <C 
ujq + uj. Then a constant energy bias Sui is added to the 
adiabatic states, and the LZ-JC Hamiltonian becomes 
H = -^^-a z - § (acr+ + aV-). The finite detuning 
Sui gives no change to the LZ process except translating 
the time when the LZ transition occurs, from t = to 
t = du/v. 

This kind of model was considered in Ref. [39], where 
highly nonclassical sub-Poissonian states were found. 
Some later work [23, 40] mentioned the LZ problem in the 
JC and Rabi models. A superconducting qubit coupled 
to a transmission-line resonator [41] can be described by 
the quantized LZ Hamiltonian in Eqs. (2) and (3) when 
the transition frequency of the charge (flux) qubit is var- 
ied linearly in time by a driving charge (magnetic flux). 
Furthermore, in this kind of circuit quantum electrody- 
namics, the strong coupling can be obtained, which al- 
lows the cases beyond the RWA [42] . 



III. LANDAU-ZENER TRANSITION 
PROCESSES, ENTANGLEMENT CREATION 
AND PHOTON DISTRIBUTION WITH THE 
RWA 

Let us assume that the initial state of the total system 
is in a direct product form, and the TLS initially starts 
from | f) and that the photon field starts from a Fock 
state, then the initial state of the whole system \ip (0)) tot 
= |t) ® \n). Under the LZ-JC Hamiltonian in Eq. (3), 
the state at time t becomes 
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A n (t) It n) + B„ (t) |; n + 1) 



(4) 



where, the time-dependent coefficients A n (t) and B n (t) 
are the solutions of a second-order Weber equation and in 
the form of combinations of parabolic cylinder functions 
[1]. Explicitly, the coefficients 



A n (t) 



± 



fJ"n,±D-l-iS n (±Z t ) . 



B n {t) = yv n ,±D_ i5n (±Z t ), 



± 
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where £>-i-i5„ {±Z t ) [D- iS „ (±Z t )] are the 
parabolic cylinder functions and the parameters 
Z t = -V2e^ /4 ^/v~/2t and S n = A 2 (n + 1) /(4v), 
because when the field mode is occupied by n photons, 
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the splitting A is enhanced by a factor y/n+ 1, as 
compared with the standard LZ model. The parame- 
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ters u n> ± and /i n ,± satisfy v n .± = =F/i„,±e ™ /4 /\AV 
Moreover, if the bias energy uio is finite, then we have 
/V+ = ^-i5„ (-^to) /[D-iS n (-Z to ) D_i^ iSn (Z to ) + 
D^is n (Z to )D-i- l5n (-Z to )} and /i„,_ = 
H ni+ D-i$ n (Z to ) /D_i$ n (—Z to ). Usually, one con- 
siders the limit cases when the bias energy vt is switched 
from a large negative value to a large positive value; 
then with the asymptotes of the parabolic cylinder 
functions, we find A n (t) ~ exp(— nS n ), which implies 
that for an initial state |f n), the final probability of the 
TLS staying in |f) is Pf >n = P-f-o exp(— irA 2 n/2v), with 
P-f-,0 = cxp(— ttA 2 /2v) denoting the vacuum case. 

When the initial state of the field is a superposition of 
Fock states, the initial state of the whole system reads 
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where the superposition parameter C„ satisfies 
Sri°=o \Cn\ 2 = 1- Then, multi level LZ transitions 
are expected at the avoided-level crossings of the 
diabatic energy levels, and we sketch the energy-level 
diagram in Fig. 1 (a). The final probability of the TLS 
staying in |t) becomes an average of Pf.n, and then the 
final LZ transition probability at t = oo reads 

P LZ (oo) = 1 - P ti0 V \C n \ 2 exp . ( 7 ) 



ra=0 
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A. LZ processes for coherent superposition states 

Hereafter we focus on the nonclassical properties asso- 
ciated with the following superposition of coherent states: 



|V(0)> 



P h 



1 



(l< 



*)), 



(8) 



where A 2 = 2(1 + costfe" 2 !"! ) is a normalization con- 
stant and, for simplicity, we assume a to be real. Such 
states are superpositions of classically distinguishable 
states and involve fundamentally nonclassical properties. 
Therefore, they are important for investigating funda- 
mental tests of quantum theory and in many quantum 
information processing tasks [32-34]. 

As shown in Refs. [43, 44], for an anharmonic oscilla- 
tor with the Hamiltonian H = v (aJa + i) +x (a^aj , an 
initial coherent state \a) will evolve into the coherent su- 
perposition states in Eq. (8) with a superposition phase 9 
corresponding to different evolution times. And at some 
special time such as t = tt/ (2%), the so-called "Yurke- 
Stoler" coherent state with 8 = tt/2 can be achieved. In 
an optical system, the coherent superposition states with 
large amplitude can be generated by using homodyne de- 
tection and photon number states as resources [34] , and 
the superposition phase 9 is related to the photon num- 
bers. 



Recalling the initial state in Eq. (6) , the superposition 
coefficient becomes 



It follows that we can calculate the final LZ transition 
probability exactly 



P LZ (oo) = 1 - 
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(10) 

It turns out that Plz (oo) now depends on the initial con- 
ditions \a\ and 9; the former is associated with the aver- 
age photon of field and the latter determines the types of 
superpositions. Obviously, for \a\ 2 = and cos 9 ^ — 1, 
then Plz (oo) = 1 — Pj-.o, corresponding to the standard 
LZ probability. Whereas in the limit \a\ 2 —> oo, with 
finite ratio A 2 /v ^ 0, Plz (oo) will tend to unity mono- 
tonically. 



1. Yurke-Stoler state 

When 9 = tt/2, the superposition state is the so- 
called Yurke-Stoler (YS) coherent state [43]: |a) Y s = 
(|a) + i |— a)) /N n n. The average photon number of 



Thus, the LZ probability for the YS state 



becomes 



P LZ (oo) = 1 - P t , exp - \a\ 2 (1 - P t , ) 



(11) 



which reveals the dependence of Plz (oo) on the ratio 



A 2 /ij and the average photon number 



Obviously, 



enhancing |a| and the ratio A 2 jv will increase the final 
LZ probability Plz (oo). 



2. Even coherent state 

When 9 = 0, the photon state is the so-called "even 
coherent state": |a), = (|a) + |— a}) /No, with iVg = 

2(1 + e _2|Q!|2 ). This state refers to the fact that the 
photon number distribution is nonzero only for even 
photon numbers with the average photon number h = 

2 \a\ 2 (l - e~ 2 l Q l 2 ) /Nq. Then the final probability 



Pi 



LZ.+ 



(oo) = 1 - P t)0 



cosh(|a;| P t , ) 
cosh lal 2 



(12) 



3. Odd coherent state 

When 9 = tt, the photon state is an "odd coher- 
ent state": |a)_ = ^ (|a) - \-a)), with A 2 = 2(1 - 
e — 2|a| ^_ £ or -p^^if.^ nly an odd number of photons have 
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Figure 1: (Color online) (a) Adiabatic energy levels of the 
quantized LZ Hamiltonian (3) within the RWA. The coupling 
A = 0.5 and frequency uj = 10 are both in units of ^Jv, also 
H = 1. The arrow marks the point where the avoided cross- 
ings are located, (b) LZ probability Plz(£) as a function of 
time, in units of l/\/v, for the coupling A = 0.5, and for vari- 
ous values of the average photon number |a| 2 . The horizontal 
black dashed lines show the analytical results in Eq. (11). (c) 
Final LZ transition probability Plz(oo) versus average pho- 
ton number \a\ 2 for various couplings A = 0.1, 0.2, and 0.5. 
(d) Final LZ probability Plz(oo) as a function of the super- 
position parameter 8, for the coupling A = 0.5 and various 
values of |q| 2 . 



a nonzero probability and the average photon number 
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/N*. Then the LZ probability 
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si nh(|a| P t , ) 
sinh |a| 2 



(13) 



When | a | approaches zero, we have Plz.- (oo) — » 1 — 
Pj? ) because the odd coherent state \a)_ tends to the 
Fock state |1). 



4- Numerical results of LZ processes for various types of 
coherent superposition states 

Let us first focus on the YS coherent state and numer- 
ically study the LZ processes. In Fig. 1 (b), we show that 
for a weak coupling A = 0.5 (in units of y/v), the LZ tran- 
sition occurs near the point t = 0, at the avoided cross- 
ings. The coherent oscillation of Pvzit) is enhanced by 
increasing the average photon number \a\ . Figure 1 (c) 
shows the final probability Plz(oo) as a function of \a\ 
for different couplings A. As expected in the analytical 
results, for increasing \a\ 2 we find that Plz(oo) tends to 



1 monotonically, and larger couplings A accelerate this 
increase. The superposition parameter 9 gives a periodic 
contribution to Plz(oo), as shown in Fig. 1 (d). For small 
average photon numbers \a\ < 2, there is a clear max- 
imum of Plz(oo) at = 7r (odd coherent state), which 
implies that not only the average photon number but 
also the proper superpositions of coherent states can en- 
hance the LZ probabilities. When \a\ is larger, however, 
there is hardly any effect of the superposition parameter 
9. Then all the superpositions of coherent states pro- 
vide the same asymptotic value of Plz(oo), as shown in 
Fig. 1(d). 



B. Entanglement between the quantum TLS and 
the photon field 

Due to the coupling terms in the LZ-JC Hamiltonian, 
the dynamics will produce entanglement between the 
quantum TLS and the photon field. The aim of this sec- 
tion is to reveal the connection between the LZ transition 
and the entanglement creation. The concept of purity 
can be employed to characterize entanglement. Based 
on the reduced density of TLSs, purity is determined 
by the linear entropy, defined by E\ (t) = 1— Trp^, (t). 
In terms of the elements of the density matrix, we have 

El (t) = 1 - Ei,j=t,I K*l Atls |j)| 2 - In the LZ process 
we find the matrix elements of the density matrix as a 
function of time 



l(t|p TLS ll>| 2 = 



l(t| Atls It) I 

IUIp-tls It) I 



[1 - Plz (t) 
Plzit). 



E 



{~\r\a\ 2n+1 A n+1 Bl 



n=0 
2 



V(^+l>! 



(14) 



Due to the exact solutions of A n (t) and B n (t), the linear 
entropy can be analytically obtained. If we choose even 
(or odd) coherent states in Eq. (8) for 9 = (or it) , there 
will be a simple form of the linear entropy 



Ei (t) - 2P LZ (t) [1 - P LZ (t)] , 



(15) 



which implies that the TLS and the photon field can 
achieve full entanglement, in the sense that, after tracing 
out the photon states, no coherence between |f) and ||) is 
left; i.e., the antidiagonal elements in the reduced density 
matrix of Eq. (14) vanish. In this case, the entanglement 
is absolutely determined by Plz (i) • At finite times, when 
-Plz (t) suddenly jumps to a nonzero value but is less than 
1/2, Ei (t) increases fast to a steady value. However, if 
Plz (t) is larger than 1/2, Ei (t) decreases. The entan- 
glement will be less perfect if other initial superposition 
parameters (9 ^ 0, tv) are chosen, because there will be 
nonzero off-diagonal elements in the reduced density ma- 
trix. The entanglement dynamics strongly depends on 
the LZ transition probability, and this is confirmed by 
the numerical results shown in Fig. 2. 
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Figure 2: (Color online) (a) Linear entropy Ei(t) as a function 
of time (in units of 1/y/v) for the coupling A = 0.5 (in units 
of y/v) and for various values of the average photon number 
\a\ 2 ■ (b) Final or long-time linear entropy Ei (oo) as a function 
of \a\ 2 , for various couplings A = 0.2, 0.5, and 1. (c) Final 
linear entropy Ei (oo) versus the superposition parameter 6, 
for the coupling A = 0.5 and various \a\ ■ 



For an initial YS coherent state, the dynamics of the 
linear entropy Ei (t) for finite times is shown in Fig. 2 (a). 
The entanglement between the TLS and the photon field 
is created when the LZ transition occurs. Analogous to 
the LZ probability, Ei (t) oscillates with time and tends 
to a steady value. However, it is not always the case that 
the larger \a\ 2 is, the stronger the entanglement will be. 
In Fig. 2 (b), we plot the final linear entropy Ei(t = oo) 
versus \a\ . Clearly, for small values of the coupling A, 
the linear entropy Ei(oo) has a maximum. However, for 
larger couplings, like A = 1, increasing \a\ only sup- 
presses -Ej(oo) rapidly. Note that in the large-coupling 
A limit, Plz (oo) achieves unity adiabatically, which im- 
plies that the final state is separable. In Fig. 2(c), the 
linear entropy £7(00) versus the initial superposition pa- 
rameter 9 for a small coupling A = 0.5 is considered. 
The odd coherent state at 9 = it has a maximal entan- 
glement for a small photon number \a\ 2 = 0.3, whereas it 
has a minimal entanglement for a larger photon number 
M 2 =2. 



C. Photon distribution of the field 

One of the best-known nonclassical effects is the gen- 
eration of sub-Poissonian (or super-Poissonian) photon 
statistics of the light field [45, 46]. A coherent state |a), 
which can be regarded as a state with the "most" classical 
behavior, yields a Poissonian distribution, i.e., the vari- 
ance of the number operator h = a^a is equal to the mean 



photon number: (An) 
the Q parameter [45], 

Q-- 



\a\ . Mandel introduced 



(An)' 



1, 



(16) 



which characterizes the departure from the Poissonian 
distribution, i.e., the nonclassical property. When Q = 0, 
the state is called Poissonian, while, for Q > the state 
is super-Poissonian. If — 1 < Q < 0, the statistics is sub- 
Poissonian. It is known that the Yurke-Stoler coherent 
state is Poissonian. However, under time evolution, the 
parameter Q changes with time as Q(t) = (An) 2 /n t — 1. 
Thus the initial Poissonian may turn into sub-Poissonian 
or super-Poissonian. It is interesting to consider how the 
photon distribution changes during the LZ process. 

For the superposition of coherent states in Eq. (8), at 
finite times, we obtain the average photon number 



n t = n + P LZ (t) . 



(17) 



where n = 2|a| 2 ( 1 — e 2 ^ 2 cos 9 J /N$ indicates the 
initial average photon number. The expression given 
above reflects that the dynamics of n t is dominated by 
the LZ transition probability Plz (t). The asymptotic 
behavior at infinite time will be 77,00 = no +Plz (oo)- We 
also obtain the average value of n 2 at t = 00, 



<n 2 > 



-A\a\'P< 



N*e\ 
\a\ 4 



3n + Plz (00) , 



(18) 



where Plz (00) is shown in Eq. (10). By the definition in 
Eq. (16), the asymptotic value of Q(t) becomes Q{oo) = 
(An) m jfiaa — 1, where the variance (An) M = (« 2 ) 00 — 
n 2 ^. We shall now concentrate on some limiting cases: 
(i) when |a| = 0, for finite LZ rate A 2 /v ^ 0, we have 
Q(oo) < 0, i.e., the LZ transition induces sub-Poissonian 
statistics in the photon field; and (ii) when \a\ 2 is very 
large, we have Q(oo) — > 0, in which case the LZ transition 
has no effect on the photon distribution. 

In Fig. 3(a), we numerically plot the Mandel param- 
eter Q(t) versus time t. Obviously, the photon statis- 
tics of the field changes suddenly when the LZ transi- 
tion occurs. Near the avoided-level-crossing point, for 
I a 1 2 > 0, super-Poissonian and sub-Poissonian statistics 
appear alternately. The final Mandel parameter Q(oo) 
versus \a\ is shown in Fig. 3 (b), for weak couplings A, 
both sub-Poissonian and super-Poissonian statistics can 
appear with increasing \a\ 2 . For the large \a\ 2 limit, the 
photon distribution finally tends to Poissonian, which is 
not shown here. We also show Q(oo) versus the LZ pa- 
rameter A/y/v in Fig. 3(c). When \a\ = 0, Q mono- 
tonically decays with A/y/v from to —1, whereas, for 
a finite average photon number such as \a\ 2 = 1, super- 
Poissonian statistics also appears. However, a large \a\ 2 
will erase the nonclassical effects revealed by the sub- 
Poissonian (or super-Poissonian). 
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Figure 3: (Color online) (a) Mandel parameter Q(t) as a func- 
tion of time in units of 1/ y/v for various values of the average 
photon number \a\ 2 and the coupling A = 0.5 in units of y/v. 
(b) Final Mandel parameter Q(oo) as a function of \a\ 2 for 
various couplings A = 0.2, 0.5, and 1. (c) Final Q(oo) versus 
coupling A for various \a\ 2 . The horizontal black dashed lines 
mark Q — and correspond to a Poissonian distribution. 



IV. LANDAU-ZENER TRANSITION WITHOUT 
THE RWA 

In this section, we consider the LZ transition without 
the RWA. A reasonable comparison between the solutions 
with and without the RWA will enable us to understand 
the contribution from the counter-rotating terms. 

A. YS coherent state case 

Without the RWA, due to the counter-rotating terms 
<7 + at and cr_a, the total number operator N = a^a+o- z /2 
is not conserved by the Hamiltonian in Eq. (2), i.e., 
[N,H] 7^ 0. Then the Hamiltonian (2) cannot be ex- 
actly diagonalized and has a fundamentally different en- 
ergy structure from the Hamiltonian (3) within the RWA, 
no matter how small the coupling A is. From the adia- 
batic energy spectrum in Fig. 4 (a) , clearly, one can find 
two groups of avoided-level crossings, which are signif- 
icantly different from the RWA case in Fig. 1 (a). The 
groups of avoided-level crossing are formed not only be- 
tween the states |f n) and \l n + 1) but also between 
|t n) and |,|,n— 1) (when n > 2), with level splittings 
Ay/n + 1 and Ay/ri, respectively. Note that the two 
groups of avoided crossings are approximately indepen- 
dent whenever the time between the successive avoided 
level crossings, i cross = 2ui/v, exceeds the duration of 
an individual LZ transition, tlz ~ max. {\ / y/v , A/v} 
[47, 48]. For our multilevel LZ problem, the couplings 




Figure 4: (Color online) (a) Adiabatic energy levels of the 
quantized LZ Hamiltonian (2) without RWA. The coupling 
A = 0.5, frequency uj = 10, and the parameters are in units 
of y/v. The vertical arrows show the points where the avoided 
crossings are located, (b) LZ probability Phz(t) as a function 
of time (in units of 1/y/v) for various values of the average 
photon number \a\ 2 , and for the coupling A = 0.5 and fre- 
quency ui = 10. The horizontal black dashed lines show the 
analytical results of Eq. (11) with the RWA which, in the mid- 
dle, approximately agree with the first-stage LZ transitions. 



become A„ = Ay/n + 1 (or Ay/n). Thus the "indepen- 
dent LZ transition approximation" holds as long as the 
Fock states \n), with n > 4uj 2 /A 2 , are not occupied, i.e., 
when uj > max{ y/v/4, \a\ A}. 

The appearance of the second set of avoided-level cross- 
ings then allows the occurrence of the second transition. 
Therefore, two stages of the LZ transition are predictable 
and numerically shown in Fig. 4(b). We choose to = 10 
and A = 0.5 (in units of y/v), and then there exist two 
(almost) independent LZ transitions. The second transi- 
tions just occur in the vicinity of the second set of avoided 
crossings at t = 2uj/v. The analytical results for Plz(oo) 
[in Eq. (11)] within the RWA are marked by dashed lines. 
For different |a| cases, the first-stage LZ probabilities 
nicely agree with the RWA results. However, after the 
second stage, the final LZ probabilities significantly de- 
viate from the RWA results. 

In order to reveal the richness of the dynamics with- 
out the RWA, in Fig. 5, we have numerically calculated 
LZ probabilities for \a\ 2 = 1 and various values of the 
couplings A and frequencies uj. Obviously, once small fre- 
quencies such as uj = 1 are chosen, the time t CIOSS = 2uj jv 
is of the order tlz , where the two LZ transitions start to 
interfere with each other. By increasing the frequency 
uj, there are two visible stages of the LZ transition. The 
first-stage LZ probabilities do not depend strongly on the 
frequency uj, as long as \a\ A <C ui is satisfied. However, 
the second-stage LZ probability is influenced strongly by 
the frequencies ui. 

If all the avoided crossings are well separated, we can 
approximately treat the transitions as being independent 
and compute the transition probabilities as joint proba- 
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Figure 5: (Color online) LZ transition probability Phz(t) as a 
function of time (in units of 1/y/v) and various values of the 
coupling A (in units of y/v) and frequency w (in units of y/v). 
The average photon number \a\ 2 = 1. The horizontal red 
dash-dot lines show the analytical results of Eq. (11) within 
the RWA. The vertical black dashed lines indicate the times 
between the successive avoided-level crossings, i.e., the two 
stages of the LZ transitions. 



bilities [28, 48, 49]. For the Hamiltonian (2) without the 
RWA, the final probability to find the TLS at |f) from 
the initial state |f n) is [49] 

P t ,„_> t = P t ,„_iP t ,„ + (1 - P t ,„_i) (1 - P t , n _ a ) , (19) 

where P t ,„ = exp[-7rA 2 (n + l)/2v] = P+q 1 , and the 
expanded form of the equation above still holds for 



the n = and 1 cases. By substituting P 



into 



Plz (oo) = 1 - X^Lol^™! p t,n->t> we obtain the LZ 
probability in the independent-transition approximation 
and weak couplings: 



Plz (oo) 



P 



t,o 



[/a,e(Pt,o) - fa,e(Pf,o)] 



(20) 



where the coefficient K a _g = e ' Q /(I + cost9e 2 ' Q ' ) de- 
pends on | a | 2 and 9. The function f a g is defined by 
fa,e( x ) = (l + x) (e^ 2x +cos9e-^ 2x y When \af -> 
0, for cosi9 ^ —1, we find that Plz (oo) tends to the 
standard LZ probability 1 — Pf,o, whereas, for cos 9 = — 1, 
Plz (oo) approaches 1 — Pjr - In the large photon number 
limit \a\ 2 -> oo, Eq. (20) gives P LZ -> 0. 

By performing a numerical time integration for £ 6 
[—50,50], in units of 1/y/v, we calculate the long-time 
LZ probability as a function of |a| for different values of 
A; these results are shown in Fig. 6 (a). Clearly, the LZ 
probability possesses a nonmonotonic behavior. When 
the coupling is sufficiently weak, such as A = 0.1, the 
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Figure 6: (Color online) (a) Without the RWA, the long- 
time LZ probability Plz as a function of \a\ 2 for different 
couplings A, in units of y/v. The frequency ui = 10 (in units 
of y/v) and the numerical time integration is performed over 
[—50, 50]. The blue solid triangles, black solid circles and red 
solid squares denote the numerical results for ui — 0.1, 0.2, 
and 0.5, respectively. The dashed lines correspond to the 
analytical results of Eq. (20), and the blue, black and red lines 
denote oj = 0.1, 0.2 and 0.5, respectively, (b) Linear entropy 
Ei{t) and Mandel parameter Q(t) as a function of time (in 
units of 1/y/v). The average photon number is \a\ 2 = 1, the 
coupling A = 0.1, and the frequency uj = 10, in units of y/v. 



LZ probability increases with |a| monotonically in quite 
a long region of |a| . However, a stronger coupling like 
A = 0.5 makes the LZ probability achieve a maximum 
value quickly. Unfortunately, the approximate result in 
Eq. (20) is not in good agreement with the numerical re- 
sults, which is unlike the thermal state case [49]. This 
is because the initial photon state is of Poissonian (sub- 
or super-) statistics, and the independent-transition con- 
dition to > max{y/v/4, \a\ A}, will be destroyed by in- 
creasing |a|. Moreover, the LZ process strongly depends 
on the frequencies to even in the weak-coupling region 
(see Fig. 5). Nevertheless, the approximation result of 
Eq. (20) still indicates a significant fact that the long- 
time LZ probability no longer increases monotonously 
with the photon number |a| . 

In Fig. 6 (b), the dynamics of entanglement and Man- 
del parameter Q can also present two transitions. From 
all the numerical results shown in Figs. 4-6, we find that 
there are some qualitative differences between the results 
within and without the RWA. Nevertheless, under certain 
conditions, the RWA results are in good agreement with 
the first-stage LZ transition. Thus one can choose prop- 
erly weak couplings A and large frequencies tu to extend 
the time interval between the two LZ transitions, when 
the RWA is valid. 

For comparison, we also consider the thermal state 
case, and below we calculate the LZ probability with the 
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is consistent with the RWA results, and the second-stage 
LZ transitions also confirm the approximate results of 
Eq. (23). Moreover, for different frequencies uj, although 
the curves differ strongly around t = 0, the LZ probabili- 
ties converge toward the same value, which is significantly 
different from the case of the YS coherent state, where 
the LZ probabilities strongly depend on w and do not 
converge to a single value (see Fig. 5). 



Figure 7: (Color online) The initial state of the photon field 
is now a thermal state. Without the RWA, we plot the LZ 
probability Plz(£) as a function of time in units of l/\/v for 
various frequencies ui = 1, 10, and 20, in units of i/v. The 
coupling A = 0.1 in units of yTj, and the scaled temperature 
T /uj = 1. The horizontal red dash-dot lines mark the analyti- 
cal results in Eq. (22) with the RWA, and the horizontal black 
dashed line corresponds to the analytical results in Eq. (23) 
without the RWA. The vertical dashed lines indicate the time 
between the successive two stages of the LZ transitions. 



RWA and without the RWA. 



B. Thermal state case 

When the photon field initially starts from a thermal 
state, the density matrix of the total system is 

Ptot (0) = |t) (tl ® \ expC-watyT), (21) 

where the partition function Z = (1 — e _ "/ T ) _1 (setting 
h,k B = 1). Within the RWA, the final LZ transition 
probability becomes 



Plz (oo) = 



I -Pi 



0.LZ 



l + nP( 



0,LZ 



(22) 



where the average photon number n = [exp(w/3) — 1] , 
and Po,lz denotes the standard LZ probability: Po.lz = 
1 — exp(— 7rA 2 /2c;). When the scaled temperature T/uj — > 
0, we find that Plz (oo) tends to the standard LZ prob- 
ability. 

Without the RWA, the approximate final LZ transition 
probability for finite temperature and weak coupling be- 
comes [49] 

Plz (oo) = ^- [/ T (P t , ) - / T (P t 2 )] , (23) 

where Gt = 1 — exp(— uj/T), and the function /t(^) = 
[1 — x cxp(— uj/T)] . For the initial thermal state, the 
probability p(n) of finding n photons is exponentially de- 
pendent on uj, which is quite different from the case of the 
YS coherent state, where the photon distribution is Pois- 
sonian and independent of uj. Consequently, in Fig. 7, for 
an initial thermal state, we plot the LZ probability versus 
time for weak coupling A. The first-stage LZ probability 



V. CONCLUSIONS 

We have investigated the LZ transition in a composite 
system of a TLS coupled to a single-mode photon field. 
The initial state of the field was chosen as a superposi- 
tion of coherent states. Within the RWA, we analytically 
obtained the LZ probability as a function of the average 
photon number and the superposition parameter. By in- 
creasing the average photon number and choosing proper 
superposition parameters, one can enhance the LZ proba- 
bility. We also found that both the creation entanglement 
(between the TLS and field) and the photon distribution 
change drastically when the LZ transitions occur, which 
is helpful for revealing the effects of the LZ transition on 
the whole quantum system. 

Beyond the RWA, we found some qualitative differ- 
ences from the RWA results. The final LZ probability 
no longer monotonically depends on the average photon 
number. In addition, two obvious stages of the LZ tran- 
sition appear in the vicinity of the successive avoided 
crossings, and the RWA results can only indicate the first- 
stage LZ probability. The final LZ probability, after the 
second stage, significantly deviates from the RWA results 
and strongly depends on the frequencies uj, even for weak 
couplings A. 

We found that the LZ dynamics is quite different from 
the thermal state case, which is due to the Poissonian 
distribution in the superposition of coherent states. Al- 
though the RWA fails in estimating the final LZ probabil- 
ity, it plays an important role in characterizing the finite- 
time coherent oscillations generated by the LZ transi- 
tions. With sufficiently weak coupling A and large fre- 
quency uj, one can extend the time interval between the 
two LZ transitions when the RWA is accurate enough. 
Very recently, the authors of Ref. [50] found the absence 
of vacuum- induced Berry phases without the RWA. Our 
results also provide examples indicating that the RWA 
leads to faulty results. 
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